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Abstract 

Recently, it was found by Greene and Levin that the Casimir energy of certain combi- 
nations of massless and massive fields in space with extra dimensions play a crucial role 
in the accelerated expansion of the late-time universe and therefore it could serve as a 
candidate for the dark energy. It also provides a mechanism in stabilizing the volume 
moduli of extra dimensions. However, the shape moduli of the extra dimensions were 
never taken into account in the previous work. We therefore study the stabilization 
mechanism for both volume and shape moduli due to the Casimir energy in M 1+3 x T 2 . 
The result of our study shows that the previously known local minimum is a saddle 
point. It is unstable to the perturbations in the direction of the shape moduli. The 
new stable local minima stabilizes all the moduli and drives the accelerating expansion 
of the universe. The cosmological dynamics both in the bulk and the radion pictures 
are derived and simulated. The equations of state for the Casimir energy in a general 
torus are derived. Shear viscosity in extra dimensions induced by the Casimir density 
in the late times is identified and calculated, it is found to be proportional to the Hub- 
ble constant. 
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1 Introduction 



According to the latest data on Type la Supernovae pQ and Cosmic Microwave Background 
Radiation (CMBR) [2j, it is strongly believed that the universe consists of a sort of vac- 
uum energy, namely dark energy, which contributes the accelerated expansion in three- 
dimensional space. Unfortunately, the exact form of the dark energy has not yet been un- 
covered until now. The prominent candidates for dark energy are the cosmological constant, 
and models of scalar fields, such as the quintessence and moduli fields. 

In the standard cosmological model where the acceleration of the universe is taken into 
account by a positive cosmological constant term, dark energy contributes largely, more 
than 70 % of the total density of the universe [2]. This number (roughly 10~ n eV 4 ) seems 
arbitrarily small and the known mechanisms, such as the popular TeV-scale supersymme- 
try (SUSY) breaking scenario or any top-down high-scale particle physics mechanisms, fail 
to produce it. 

In recent years, theories with large extra dimensions have received an explosion of interests 
as they provide new solution to the hierarchy problem. Recently, it was found that Casimir 
energy of massless and massive fields embedded in higher- dimensional spacetime could play 
a crucial role of dark energy with additional significant properties [31 H]. The Casimir energy 
not only drives the expansion of universe acceleratedly, but also stabilizes the volume moduli 
of extra dimensions. However, the shape moduli, T\,T2, were not included in the work 
of Greene and Levin. In this work we therefore take into account these moduli in the 
cosmological dynamics by assuming that the extra dimensions are T 2 . The phenomenological 
implications of nontrivial shape moduli were pointed out in [5j El [7] . Shape moduli can have 
dramatic effects on the Kaluza-Klein spectrum, for example, they can induce level-crossings 
and varying mass gaps. They can also help to eliminate light KK states. It should be 
interesting to investigate the role of shape moduli in cosmology. 

Our work employed the calculation of Casimir energy in the non-trivial space M 4 x T 2 . 
The Casimir energy is the vacuum energy contributed from the quantum fluctuation of fields 
which satisfy certain boundary conditions. In fact, the Casimir energy in various spaces 
including a distorted torus was studied in earlier works [HI [9l [10l H] . The standard approach 
for determining the Casimir energy is the zeta function regularization [TT] . 

Our result shows that the minimum of potential in the previous work [3 J (ti = 0, r 2 = 1) 
was the unstable local minimum while the true local minimum locates at specific points in 
the moduli space, r\ = ±1/2, r 2 = v3/2, confirming the result of Ref. |@|. At this local 
minimum the potential stabilizes all moduli and also sources the accelerated expansion of 
the four dimensional universe. 

This paper is organized as follows. In Section 2 we review cosmological dynamics on 
M 1+n x T p spacetime. In Section 3 we present the mathematical calculation to determine the 
Casimir energy of massive and massless fields in the spacetime with toroidally compactified 
extra dimensions. Then we go on to construct effective potential contributed by Casimir 
energy of massive and massless field in M 1+3 x T 2 spacetime in section 4. The numerical 
evidences of the stability of moduli space are presented in section 5. In section 6 we present 
our conclusions. 
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2 Cosmological Dynamics in M 1+n x T p 

Our study of cosmological dynamics is based upon the application of Einstein's general 
relativity on the product space M 1+n x T p , between a (1 + n) -dimensional spacetime and 
a p- dimensional toroidally-compactified space. As a whole, the total number of spatial 
dimensions is d = n + p. We assume the cosmological ansatz 



ds 2 = g iiV [x)dx il dx u + hij(x)dy' l dy\ 



where the metric represent the p-dimensional compact space with i,j = 1, • • • ,p and g^ u 
for the (1 + n) -dimensional noncompact spacetime with p,, v — 0, • • • , n. Let's assume also 
that the metric only depends on the noncompact coordinates x^ . The compact coordinates 
are < y i < 2vr. 

In this paper, we focus our effort on the cosmological dynamics of a 4-dimensional space- 
time with two extra dimensions (n — 3 and p — 2). The metric of two-dimensional torus T 2 
takes the form 



(hij) 



T 2 



1 



Ti 



(2) 



where r = T\ + ir 2 is the complex structure (or shape moduli) and b 2 is the Kahler structure 
(or volume moduli). In cosmology, it is customary to write g^ u = a 2 (t)r]^ and = hij(t). 
In the next sections, we will assume that Casimir energy in compact direction, p(d+i)D, plays 
the roles of the dominant energy content in the universe. By using Einstein equations in 
(1 + 5)-dimensional spacetime, we obtain the following equations governing the cosmological 
dynamics: 



H a + 3H 2 + 2H a H b 



H b + 2H 2 + 3H a H b 



_ 8ttG 
8%G 



2P6D + 



mi + H 2 b + 6H a H b - —s(ft + f 2 ) = 8vrGp 6D , (3) 



At, 



1 - 



2r 2 



bd b p 6D - 2r 1 d T1 p &D H d T2 p 6D } , (4) 

2t 2 

bd b p 6D - 2r 1 d Tl p 6D H -d T2 p &D } , (5) 

r 2 



r, ( 3H a + 2H b - 2^ \ t x = -16nGr 2 |^<9 fc p 6D + 2d T1 p 6D - ^-d T2 p eD \ , (6) 



T 2 . T l - T 2 



+ 



T2 



br 2 



r 2 



+ 3H a — + 2H b - = 8nG { -^d b p 6D + 2T 1 d n p 6D ~ 2r 2 



r 2 



T2 



T2 P6D > • (7) 



where G is the 6-D gravitational constant. We have defined the Hubble constants H a = a/a 
and H b = b/b, where a dotted quantity represents the corresponding time derivative and p 6D 
is the casimir energy density in six dimensional spacetime. 
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2.1 Dynamics in the Radion picture 



Equations of motion ©-([7]) can be obtained by varying the ci+l-dimensional Einstein-Hilbert 
action: 



S= I d 1+n xd p y^-gh\^—n {1+d) - p {1+d)D (W) 



M: 



d-l 



16tt 



(8) 



with n = 3 and p — 2, where P(i+d)D(h 1 ^), H(i+d) and M* are the Casimir energy density, 
Ricci scalar and the Planck mass in (1 + d) -dimensional spacetime respectively. For later 
purpose, it is useful to perform KK-dimensional reduction of the above action from (1 + d) 
to (1 + n)-dimensional spacetime and Weyl rescaling g^ UB = Q : ^g fll/ ; Q = M^Vp/m^ 1 , 
the action takes the form 



S 



n-1 



pl 1 "* v( 1 - V u In v^V, In + \v u h ij V v h 



n 



(9) 



Note that the subscript E denotes the Einstein frame variables. Here, V p — J d p yVh = 
{2ixb) p = l p is the (invariant) volume of extra dimensions, m p i and U(h' lJ ) = VL T ^V p p(i + ct )D {h %:) ) 

1 + n . . 

Vt 1 -™ P(i +n )D\h i:i ) are the Planck mass and the effective potential in 1 + n-dimensional space- 
time respectively. We can also take p{i +n )D{h^) to be the Casimir energy density in (1 +n)- 
dimensional spacetime. 

Since we are interested in the n = 3, p = 2 case, by using the metric of two-dimensional 
torus defined in Eqn. (T5]), the action in Eqn. (J9j) can be written as 



S = d x v / -g E 



m 



P l 



16n 



n E --g E u (V^V u ^+e- 2 ^V^ 1 V u <p 1 +V^ 2 V u <p 2 )\ -U($, fa, 2 ; 



(10) 

where ip = 2\/21n&, 0i = Ti, and 02 = lnr2. Such action gives rise to the following set of 
equations: 

6iJi-i(^ + e- 2 ^ 1 2 + 2 2 ) - 



i> + 3H E ^ 
0i + 3if s 0i - 20102 
+ 3^02 + e- 2 ^0! 2 



and 



16 » 






(11) 


m pl 




16vr dU 


(12) 


m 2 pl di\) ' 


16^ dU 


(13) 


m 2 pl " 90i ' 


16vr dU 


(14) 




m 2 pl 902 ' 




)=0. 


(15) 



4H E + (^ 2 + e-^0! + 2 " 
Note that H E = (da E /dt E )/a E is the Hubble constant in the Einstein's frame. 
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3 Casimir energy in M 1+n x T p 

In this section, we will undergo the mathematical formulation to determine the Casimir 
energy, E cas , associated with a scalar field of mass M in a M 1+rt x T p space. The fermionic 
degree of freedom will contribute to the Casimir energy with the same expression except 
for an extra minus sign. We then focus on the result from our phenomenological study 
(n = 3,p = 2). 



3.1 Casimir-Energy Calculation 

Let V n = L n be the spatial volume of non-compact spacetime, and V p = V be the volume of 
compact space. If we assume L>I, the zero-point energy of scalar fields in M 1+ " x T p can 
be evaluated by 

■~ 1 L r +co i 

E cas = -(—)»£/ d n k^5" b k a k b + hv ni n 3 +M2, (16) 

where k a \ a = l,...,n is the momentum in each non-compact spatial direction, n, G Z; 
i — 1, ... ,p is the momentum number in each compact direction. 

Using the property of integration in Appendix |A] and changing variable of integration as 
v = k 2 1 '{h^riirij + M 2 ), we can express the Casimir energy as 

E C as = ^T ¥r i^Y,^ n ^ + M2 )^ dw^VTT^. (17) 

2 2tt T[n/2) f-' J 

We can convert the integral into the Gamma function by using the formulae in Appendix |A| 
as a consequence, we obtain the Casimir energy in a simple form 



2 l' Tr^rr-i) 



In our case, the compact space is T 2 and is the inverse metric from Eqn. (jSJ); therefore, 
our next task is to regularize the infinite summation in the Eqn. (jT8l) 

which is known as extended Chowla-Selberg zeta function [9j. It is worth noting that V^, = 
Z 2 = (27r6) 2 in this case. 

After a few steps of analytic manipulation by using Poisson resummation and property 
of the modified Bessel function, we obtain 



i^;^,-|^,7^;M 2 ) = b^{2T^ EH {s,r 2 b 2 M 2 ) + 2^%-^rl~X EH {s -1/2;^) (20) 
b 2 r 2 b l r 2 b 2 r 2 T(s) r 2 



^ 8tt s _ i cos(27rr 1 mfc) r ^ . / „ 6 2 M 2 
+ > -^\pT2k 2 — , - — if, i (2-KT 2 k\ m 2 H )}, 



m,k=l 
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where the Epstein-Hurwitz zeta function (eh(s; q) is expressed as 
(EH{s;q) = ^V + g)- s 



2- 

nGZ 

v ) n=l y ' 

where the prime at the first sum indicates that the term n = is excluded. A similar 
expression which manifests the periodicity of the Casimir energy with respect to T\ is also 
given in Ref. |12j . 

The expression serves as an analytic continuation of the Casimir energy where s is ex- 
tended from positive to negative values. Inserting Eqn. (1201) into Eqn. ffT8l) and eliminating 
the infinite terms due to the pole of T(s = —2) and T(s — 1 = —3) in this case, we con- 
veniently reached the regularized Casimir energy. The dropped divergent terms correspond 
to the constant total energy and the constant energy density in the bulk. Both of them do 
not depend on any parameters of the torus and therefore can be safely eliminated from the 
physically relevant Casimir effects by renormalization. The final regulated Casimir energy 
density pih 1 ^) in (1 + 3) -dimensional spacetime can then be expressed as 



P4D(b 2 ,T 1 ,T 2 ) 



Er. 



Vm 

oo 

-(47r 2 6 2 ) s |2r 2 s (r 2 6 2 M 2 )-f + ^F-^ s _i(27rA:6M v ^) 



fc=i 

oo 



^~ 2 k=l 



a i— 1 5-i cos(27rriA;m) / b 2 M 

> k 2 — . v ; — K s _i (2nkT 2 \ m 2 + 

^ l m 2 + miy-k 2 V r 2 



k,m=l 



(22) 



In the case of massless scalar fields (M = 0), the Casimir energy density becomes 

P,n(b 2 , n, r 2 ) = -{A^b 2 ) s {r^Y{^ - S )C(1 - 2a) + ^^"^(l - s)((2 - 2s) 

oo ^ 

+4^ (— ) S "^cos(27rmA;ri)K s _i(27rmA;r2)}. (23) 



The Casimir density in (1+3+2) dimensions is given by p%r> = p^D / (2nb) 2 . 

As it is pointed out in the work of Ponton and Poppitz |4J. Since the symmetry r — > 
— 1/t, r^r + lof the torus is preserved in the Casimir energy expression, it is sufficient 
to consider only the fundamental region where r > 1,-1/2 < t\ < 1/2 of the shape 
moduli space. In the fundamental region, there are two minima and one saddle point of the 
magnitude \p\ of the Casimir energy density. The saddle point locates at r± — 0, r 2 = 1 and 
the two minima locate at T\ = ±1/2, r 2 = VS/2. This is shown in Figure [lj 
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3.2 Analysis for small bM 

In the limit of bM C 1, we recalculate the Casimir energy by performing the binomial 
expansion with respect to small bM before regularization, and keep only the leading-order 
terms. It can be demonstrated that the process of regularizing each term after performing 
binomial expansion is NOT equivalent to the process of regularizing the whole expression at 
once if s = —2 is set beforehand. When we set s = (1 — d)/2 — —2, the binomial expansion 
of Eqn. (fl9|) gives only three terms with orders of (bM)°, (bM) 2 , and (bM) 4 , whereas the 
regularization of the full expression before setting s = —2 as in Eqn. (120]) , which gives 
Eqn. (|22|) as a result, generically leads to an infinite series of bM, even after setting s = —2 
in the final expression. 

Without setting s = — 2 before regularization, the precise dependence of the coefficients 
of the 6M-binomial expansion to the moduli parameters r 1; r 2 will be determined. The small 
bM expansion is obtained subsequently. 

We begin by replacing with the form of the inverse metric of T 2 in Eqn. ffl8l) and 
using Mellin transform (see Appendix [AJ 

cas " 2 ( L } ^r(-i)„f c Jo 



2' ni,n,2SZ 



1,2tt, 1+2s T(s) 



2 y L' Tx^Yi , 

\ 2> m,n2& 



'> J 1—11 f-l-i TtnCZ.//. s 



where the second line is obtained by changing the dummy variable v = t/b , and the final 
line is obtained by expanding the Taylor series for e~( bA/ ) 2 . We can determine the double 
summation in Eqn. (|2~4"I) by using the result in Eqn. (TEH) . (|2"UI) ; as a consequence, the Casimir 
energy density in five spatial dimensions takes the form, 



(24) 



( 1 \j h 

/06o(6 2 ,Ti,r 2 ) = -(47r 2 6 2 ) s - 1 ^^ T ^(6M)^{47r J V^ Yl (-) S+i ^ ^irmk^K^i (2nmkr 2 ) 

j=0 3' m.k=\ m 

+ ^-hr s 2 ^T(\ - a - j)C(l -2s- 2j) + ^-Vf— ^(1 - a - j)((2 -2s- 2j)). 
In the limit bM -C 1 for s = —2, the Casimir energy density then becomes 

p eD (b 2 , 7i, r 2 ) ~ - ^ {d - C 2 (6M) 2 + C 3 (6M) 4 } (25) 
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where 



Ci = 



7T 



f r 2 " 2 r (I) C(5) + vr- 3 r^r(3)C(6) + 4^ £ 





cos (2irmkTi ) K_ 5 / 2 (27rmfcT 2 ) 



C 2 = 




C(3) + 7r- 1 r 2 2 r(2)C(4) + 47r v ^ £ 




cos(27rmfcr 1 )ii'_3/2(27rm/cr2), 



C 3 = 



£r 2 r(l)C(2) + 2ttV* £ (^) 



cos (27rm/cri ) !/ 2 (27rmfcT2 ) . 



m,fc=l 



(26) 



In the next section, the total Casimir density for small bM and the full expression will be 
numerically compared. The true minimum of the potential, induced from the Casimir energy 
density located at a point (n, r 2 ) = (±1/2, y/3/2), appears only when the full expression is 
evaluated. 

4 Particle spectrum and effective potential for moduli 
fields 

It is demonstrated in Ref. [I] and Ref. [3] that a careful mixing of massless and massive, 
bosonic and fermionic degrees of freedom of the bulk fields can lead to a Casimir energy den- 
sity with local minimum with respect to the scale factor, b, of the compact extra dimensions. 
In the torus case with the shape moduli ti,t 2 , it can be shown that the true minimum of 
the mixed Casimir energy density (and thus the potential) locates at T\ = ±1/2, r 2 = v / 3/2, 
in contrast to the case of undistorted torus considered in the previous work where the shape 
moduli are set to T\ = 0, r% — 1. 

The simplest model of the bulk fields in our M 1+3 x T 2 space consists of a massless boson, 
a massless fermion, a massive fermion with mass M, and a massive boson with mass AM. 
It was found that for the range 0.40 < A < 0.42 and M = 5, the mixed Casimir density 
has local minimum with respect to the scale factor b, and the moduli T\, t 2 . Since the mass 
of the boson is different from the mass of the fermion, this is the scenario where SUSY is 
broken in the bulk if it exists at higher scales. There is no particular reason for why the 
ratio of the masses of the massive boson and fermion took the specific value in this range. 
If it has anything to do with SUSY breaking, it is desirable that we are able to establish 
a SUSY breaking mechanism where this specific ratio of the masses A could be explained 
or distinctively selected. From phenomenological point of view, it is desirable that these 
massless and small-mass bulk fields are sterile neutrinos for they can explain the smallness 
of neutrino masses in four dimensions. For further details, see Ref. [i~3].[14j. 

An important issue in mixing bosonic and fermionic degrees of freedom to obtain the 
total Casimir energy density with a local minimum is the positivity of the energy density. 
Generally, the value of the total Casimir density at T\ = ±1/2, r 2 = V^/2 is lower than the 
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value at the saddle point T\ = 0, r 2 = 1, for all range of A. However, for certain ranges of 
^ i e -9- ^ 0.407), the density becomes negative around the true minimum and therefore 
violates the positive energy condition. A negative value of the density will not stabilize 
the dynamics and the size of the torus. We therefore choose the value A = 0.408 for our 
simulation of the cosmological dynamics. Figure [2] shows the total Casimir energy density 
for the spectrum of massless and massive particles mentioned above. 

The plot of the total Casimir density in (l+3+2)-dimensional spacetime using the full 
expression, Eqn. (1221) . in comparison to the plot from the small bM approximation, Eqn. ([25]) . 
is given in Figure [31 The true minimum at T\ = ±1/2, Ta = Vo/2 only exist in the full 
expression case. This can be understood considering b m i n M rs 0.67 and is somewhat close 
to 1, resulting in a bad approximation of the expression due to higher powers of bM being 
neglected. It is therefore required that we use the full expression of the total Casimir energy 
density in the simulation of the cosmological dynamics. 

5 Evidence of stability of the moduli space and cosmo- 
logical dynamics 

By numerically solving the field equations in section [21 the stabilization of the torus and the 
accelerated expansion of large 4-dimensional spacetime can be demonstrated to occur at the 
true minimum of the Casimir energy density in the moduli space. The point t\ = 0, r 2 = 1 
is a saddle point and it is an unstable equilibrium of the dynamics. 

The rolling of the universe to the true minimum of the Casimir density is illustrated in 
Figure [4}[7J When the cosmological dynamics is initiated even within a small vicinity of the 
saddle point, T\ = 0, r 2 = 1, of the Casimir energy density, it will roll down to the true 
minimum at n = ±1/2, r 2 = x/3/2 even with minimal amount of perturbations. This is 
shown in Figure HI El Observe that it tends to roll along the trail r = 1 in the moduli space. 

When the tossing initial conditions are at a distant away from the saddle point and the 
true minimum, certain sets of the initial conditions still result in the stabilization of the 
torus moduli, ti,t 2 , and the scale factor, b, of the extra dimension as is shown in Figure El 
[71 Naturally, as long as the Casimir energy density at the stabilized value is positive, the 
acceleration of the scale factor, a, of the 4-dimensional spacetime is guaranteed. The positive 
Casimir density serves as the positive cosmological constant. 

A natural consequence of the Casimir energy that is independent of the scale factor, a(t), 
of the large dimension is the fact that it leads to w a = — 1 for the pressure p a = w a p. For the 
pressure in the compact extra dimensions, we can start by considering pb = —d(pVb) / dVf, = 
WbP, Wb of our Casimir energy density is then given by 

where p is the total Casimir energy density. Due to the dynamics of shape moduli (or Casimir 
"viscosity" in the compact space, see Appendix [Bj) , the value of Wb at the stabilized radius 
at the true minimum is fractionally smaller than —2 (around —2.16) as is shown in Figure [5] 
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A more appropriate definition of physical pressures in the distorted torus is 



Pk = T«, (28) 
where K = 4,5. This definition gives the following expressions for wk = Pk/ Pi 

h r 2 - t 2 2ri 1 r 2 - r 2 

^4 = + + )^ (29) 

2p r 2 z p p r 2 

"6 = + (30) 
^p r 2 p r 2 

By directly solving the equations of motion in six dimensions at the stabilized point where 
H a — _H& — Hf, — j\ — fi — 0, it can be shown that = —2, as is confirmed numerically 
in Figure El It is interesting to note that the value of becomes —2 at both the saddle 
point and the true minimum where the dynamics is stabilized. 

The difference of the two definitions of pressure originates from the shear viscosity induced 
by the Casimir energy in the off-diagonal components of the stress tensor. From the equations 
of motion of the 6-D universe with viscosities, Eqn. (152]) in Appendix [Bl shear viscosity at 
the stabilized point r]f l tab can be identified to be 

^stab _ 3H asta b 

P6D,min (Qr>\ 
a, stab 

where H ajSta b is the Hubble constant of the expanding four dimensions at the stabilized point 
of the compactified space. Note that we can evaluate Eqn. (Ej), ( 1501 and (15T1) at the stabilized 
point and use the definition of r]l tab to analytically confirm the numerical results in which 
w 4j5 = —2 at the stabilized point. 

We should mention here that the time scale, t s , of the simulated figures is given by 

t. = ^-H, (33) 

where b s is the scale of b, and b min ~ 0.13286 s as a result of numerical simulation. If we 
require that the stabilization time ~ 10i s is less than the age of the universe, 10 10 years, this 
will put constraint on the size b m i n of the extra dimensions T 2 , 



bmin < 0.7pm. (34) 

This is about few hundred times stronger than the constraints from table-top experiments |15j . 

It is interesting that in this kind of cosmological model, the constancy of the 4-dimensional 
gravitational constant, G4 = G/4n 2 b 2 = 1/m 2 ^ up to the early times of the universe will 
give a very strong constraint on the size of the compactified extra dimensions. Any future 
observations of the universe from very early epoch could possibly put constraints on the 
inconstancy of the gravitational constant. Such constraints will put very strong limits on 
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the size of compact extra dimensions in this kind of model where oscillatory behaviour is 
significant in the early times. 

Another important aspect of this model is the relationship between the effective cos- 
mo logical constant in 4-dimensional spacetime, A4 = 8irG4p4D,min, and the size of extra 
dimension, b min , 

A 4 = SnG^PiD^in (35) 

= 3#i )Sta6 . (36) 

This leads to the typical value of b min ps 2.4 jttm for p vac ps 1CT 11 eV 4 . The value of the 
effective size of extra dimensions, 2nb m i n ps 15/im, yields the quantum gravity scale in the 
bulk, M* ps 12 TeV. 



6 Conclusions and discussion 

The stabilization of compact extra dimensions and the acceleration of the other 4-dimensional 
part of the spacetime can be simultaneously described by the dynamics of the Einstein field 
equations in the bulk spacetime. The acceleration of the 4-dimensional "universe" occurs 
naturally once the scale of the compact dimensions is stabilized and the density of the 
Casimir energy in the bulk becomes a (positive) constant at that stabilized value. As a 
result, the apparent positive "cosmological constant" that we seem to observe in the four 
dimensional visible universe is effectively induced. This is demonstrated beautifully in the 
work by Greene and Levin [3] when the Casimir density of the undistorted torus satisfies 
w a = —1, Wb = —2 condition. 

Shape moduli of the torus can be added to the model. The true minimum of the 
Casimir energy density of the torus with shape moduli is demonstrated to be located at 
T\ = ±1/2, t% = Vo/2. The cosmological dynamic shows that a minimally small perturba- 
tion to the saddle point rolls the universe down to the true minimum. Other initial conditions 
also suggest that the universe tends to roll around r = 1 contour to reach the true minimum. 
Note that it is also possible to stabilize the moduli at the saddle point t\ = 0, r 2 = 1 but 
the initial conditions of the shape moduli fields must be fine-tuned so that t\ — 0, j\ — 0. 
Some extra-mechanisms such as Brandenberger-Vafa mechanism in string gas cosmology 
[T6] is needed for this purpose. However, as it was pointed out in [17], the stabilized point 
T\ = ±1/2 and r 2 = v^3/2 is also the fixed point of T-duality and the the enhance symmetry 
point hence Brandenberger-Vafa mechanism could also set the initial value of the moduli 
precisely to be at the stabilized point. 

The shear viscosity in the extra dimension is determined to be proportional to the Hubble 
constant at the stabilized point, T]b = 3H a ^ ta b/16iiG. Through the Einstein field equations, 
this Hubble constant of the 4-D universe is determined by the value of the Casimir energy 
density at the stabilized point. The effective four dimensional cosmological constant is also 
given by 8nGp 6Djmin . 

In this kind of model, there is a relationship between the size of the compact dimensions 
and the observed four dimensional cosmological constant. This remarkable connection is 
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induced by the nature of Casimir energy density which depends on the size of the compact 
dimension, resulting in A 4 ~ b^ in . 

It is equally important to note that the constancy of the 4-D gravitational constant up to 
very early time of the universe will provide strong constraint on the size of extra dimension in 
this particular cosmological model which expresses oscillatory behaviour at the early times. 
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A Useful Formulae 

Phase Space Integration 



j d d kf(k) 



J k^fitydk. (37) 



r(d/2) 
Poisson Resummation 

J2f(n) = f(k) = y/to^ffr™), (38) 



where 

i r+oo 

1 / „, , A kx 



f(k) = -= j f(x)e^ x dx. (39) 



If f( x ) = e - a (*+ c ) 2 , then f(k) = -^^le~ ! t +ikc . 
Integral representation of Gamma function 

/•+oo 

T(x) = / e~H x - l dt. (40) 
Jo 

The integral representation of the modified Bessel function of the second kind 

where \arg(z)\ < f ,Re(z 2 ) > 0. 
Mellin transform 



1 f°° 

= — — / dt e-'H 3 - 1 ; Re(z) > 0, Re(s) > 0. (42) 
r(s) Jo 
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B Energy momentum tensor of viscous fluid 



Let U A = (1, 0, 0, 0, 0, 0) be the 6-velocity of the cosmic fluid in comovine coordinates. In 
terms of the projection tensor Kab = 9ab + UaUb, the general energy momentum tensor of 
fluid with bulk viscosity ( and shear viscosity rj is given by: 

T AB = pU A U B + (p - (0)h AB - 2 Wab . (43) 

Here 9 = V aU a is the scalar expansion and oab = b^h^V \cUd) ~ \hAB^ is the shear tensor. 
By using metric defined in Eqn. (TjQ) and (j2J), we can show that 

= ~P, (44) 
T 1 = Ti = Tl= Pa (45) 

Tt = (p b -C b 0)-2r]AH b -H a )-^--^) (46) 



-p, 




rp2 

1 2 - 




(Pb 


- C b o) - 


(Pb 


- ( b e) - 


M 


nf 2 ^ 






r 2 




h 


~Vb 






T 2 
T 2 



-5 y ' 2r 2 2r 2 

r 3 /rr 7-7- \ , t~2 , Tin 

r-2 



T 5 5 = (p b -( b e)-2 Vb [^(H b -H a ) + ^- + ^} (17) 



Tl = 2 % [^ + (r 1 2 -r 2 2 )^|] (48) 
Tl = -Vb~ 2 (49) 

T 2 

Here we assume there is no viscosity in noncompact large dimensions (( a = rj a = 0). Ein- 
stein's equations, Eqn. (jl])-©, can be written in terms of bulk and shear viscosity as 



H a + 3H 2 a + 2H a H b = ^ |p 6D + Pa - 2{ Pb - ( b 6) + - H a ) 



H b + 2H 2 b +3H a H b = ^\p 6D -3p a + 2(p b -C b 9)-^r, b (H b -H a ] 



(50) 



(51) 

h + \3H a + 2H b - 2^J n = 16vrG {r] b h} , (52) 
fl + ^l + 3Ha ^ + 2H b ^ = 48^(^1. (53) 

T 2 T% T 2 T 2 I T 2 J 

The conservation of energy is 

p 6D + 3H a (p 6D + pa) + 2H b (p 6D + p b ) 

+ &b-6(b)H a H b ^ r]b + 4( b )H b 2 -r ]b (^ + ^) = 0. (54) 

D Try Try 
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x 2 



Figure 1: The magnitude of the Casimir energy density, \p4d\, in four dimension per degree 
of freedom for M = 5, b = 0.133. 




b 



Figure 2: The total Casimir energy density in six dimension for mixture of massless and 
massive fields for M = 5, A = 0.408, and r = a/ r\ + r| is fixed to 1. 
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Figure 3: The Casimir energy density in six dimension from small bM approximation in the 
upper figure in comparison to the full expression in the lower figure. Both are evaluated at 
their corresponding b m i n . 



16 




17 




18 



Figure 6: Rolling dynamics from other initial condition I. 




